CFD simulation of airflow over a regular array of cubes. Part II: analysis of spatial average properties     

Alberto Martilli  Jose Luis Santiago 《Boundary-Layer Meteorology》2007,122(3):635-654

In the first part of this study, results of a computational fluid dynamics simulation over an array of cubes have been validated against a set of wind-tunnel measurements. In Part II, such numerical results are used to investigate spatially-averaged properties of the flow and passive tracer dispersion that are of interest for high resolution urban mesoscale modelling (e.g. non resolved obstacle approaches). The results show that vertical profiles of mean horizontal wind are linear within the canopy and logarithmic above. The drag coefficient, derived from the numerical results using the classical formula for the drag force, is height dependent (it decreases with height). However, a modification of the formula is proposed (accounting for subgrid velocity scales) that makes the drag coefficient constant with height. Results also show that the dispersive fluxes are similar in magnitude to the turbulent fluxes, and that they play a very important role within the canopy. Vertical profiles of turbulent length scales (to be used in k–l closure schemes, where k is the turbulent kinetic energy and l a turbulent length scale) are also derived. Finally the distribution of the values around the mean over the reference volumes are analysed for wind and tracer concentrations.  相似文献   

北京市朝阳区沙子营垃圾污染物在潜水含水层中迁移规律的现场模拟实验     

刘长礼詹黔花  张凤娥张云 张胜叶浩  侯宏冰张明裴丽欣  姜建梅 《地质通报》2007,26(3):305-311

选择北京平原区水文地质、环境地质等方面都比较典型的沙子营垃圾堆放场,建立了其水文地质模型。在充分收集资料、分析得出部分计算所需参数后,采用现场弥散实验、勘查取样测试等方法,求得了该含水层的弥散系数等参数;用二维非稳定流溶质运移方程对污染物在此含水层中的迁移扩散规律、速度和污染范围等进行了模拟计算;采用现场钻探、取样测试分析等方法,评价了该含水层的实际污染状况。实验模拟计算和现场调查结果表明:污染物在潜水含水层中的运移规律遵循二维非稳定流场中的溶质运移方程,污染物在潜水含水层中的运移速度约为86.25m/a,迁移扩散主要发生在地下水流向上,侧向扩散宽度极小,是地下水流向上的1/17。  相似文献   

Migration and dispersion of trace elements in the rock–soil–plant system in areas underlain by black shales and slates of the Okchon Zone, Korea     

Jin-Soo Lee  Hyo-Taek Chon  Kyoung-Woong Kim 《Journal of Geochemical Exploration》1998,65(1):61-78

The Okchon black shale in Korea provides a typical example of natural geological materials enriched with potentially toxic elements. The Chung-Joo, Duk-Pyung, Geum-Kwan, I-Won, Bo-Eun and Chu-Bu areas are underlain by these black shales and slates of the Guryongsan Formation or the Changri Formation, which are parts of the Okchon Group in the central part of the southern Korean Peninsula. In order to investigate the enrichment levels and dispersion patterns of potentially toxic elements in the rock–soil–plant system, environmental geochemical surveys were undertaken in the above six study areas in the Okchon Zone. After appropriate preparation, rock and soil samples were analyzed for potentially toxic elements by instrumental neutron activation analysis (INAA) and inductively coupled plasma atomic emission spectrometry (ICP-AES), and plant samples by atomic absorption spectrometry (AAS). In particular, Ba, Cd, Mo, V and U in Okchon black shales are highly enriched, and their mean concentrations are significantly higher than those in black slates. These elements are geochemically associated, and might be enriched simultaneously. The highest mean concentrations of 42.0 μg g⁻¹ As, 2100 μg g⁻¹ Ba, 10.9 μg g⁻¹ Cd, 213 μg g⁻¹ Mo, 83 μg g⁻¹ U, 938 μg g⁻¹ V and 394 μg g⁻¹ Zn are found in black shales from the Duk-Pyung area. Mean concentrations of As, Mo and U in soils overlying black shales occurring in the Duk-Pyung area (30 μg g⁻¹ As, 24 μg g⁻¹ Mo and 50 μg g⁻¹ U) and Chu-Bu area (39 μg g⁻¹ As, 15 μg g⁻¹ Mo and 27 μg g⁻¹ U) are higher than the permissible level. Enrichment index values of the six study areas decrease in the order of Duk-Pyung > Chu-Bu > Bo-Eun > Chung-Joo > Geum-Kwan = I-Won areas. Relationships between trace element concentrations in soils and plants are significantly correlated, and the biological absorption coefficients (BAC) in plants are in the order of Cd > Zn = Cu > Pb, which suggests that Cd is more bioavailable to plants than the other elements. Cadmium concentrations in plant species decrease in the order of chinese cabbage > red pepper > soybean = sesame > rice stalk > corn > rice grain. From the result of sequential extraction analysis of soils, relatively high proportions of Cu, Pb and Zn are present as residual fractions, and that of Cd as non-residual fractions. Cadmium occurs predominantly as exchangeable/water-acid soluble phase in soils, and this is in agreement with the findings of high Cd concentrations in plants.  相似文献   

Modelling dispersion in complex open channel flows: Fuzzy calibration (2)     

B. G. Hankin  K. J. Beven 《Stochastic Environmental Research and Risk Assessment (SERRA)》1998,12(6):397-412

Part 1 of this study assessed the effectiveness of dispersion models based around random particle tracking (RPT) applied an overbank flow in the Flood Channel Facility at HR Wallingford. Several of the RPT models that were used predicted the mixing behaviour of the tracer in the complex channel flow to within an acceptable accuracy. If there is no uniqueness in the model structure which can be used to represent the true system to within the limitations of the available observations, then this implies an inherent degree of ambiguity in our knowledge of the physically based model structure. This suggests that we should be less forthright in the optimisation of each individual model structure, and perhaps investigate more of the parameter combinations for each model which yield feasible simulations of the system. An alternative fuzzy calibration technique is introduced which avoids the optimisation process and takes account of uncertainties in the model structure, parameter sets and observed data in prediction.  相似文献   

Dynamics of upper tropospheric stationary wave anomalies induced by ENSO during the northern summer: A GCM study     

R. Krishnan  C. Venkatesan  R. N. Keshavamurty 《Journal of Earth System Science》1998,107(1):65-90

Ensemble seasonal integrations are carried out with the COLA GCM, with a view to understand the dynamical connection between warm SST anomalies in the equatorial central-eastern Pacific Ocean and the upper level stationary wave anomalies seen during drought years over the Indian summer monsoon region. In addition, experiments with and without orography are performed in order to examine the role of the Himalayas in modulating the El Niño induced stationary wave anomalies over the summer monsoon region. The GCM simulations show a statistically significant weakening of the summer monsoon activity over India in response to the SST forcing in the equatorial Pacific Ocean. This weakening of the summer monsoon appears to be largely related to modifications of the local Hadley and Walker cells over the summer monsoon region. In addition, it is seen that the anomalous ENSO divergent forcing over the tropical Pacific Ocean can act as a potential source for Rossby wave dispersion. Here one finds the possibility of meridionally propagating Rossby waves, which emanate from the ENSO forcing region, to interact with the subtropical westerlies and generate anomalous highs and lows in the subtropics and extratropics. The quasi-stationary perturbations seen over west Asia, Pakistan and northwest India during drought years, seem to be generated by the above mechanism. An alternate mechanism that could be important for the persistence of the quasi-stationary perturbations seems to be based on the dynamic excitation of middle latitude normal modes which can extract energy from the zonally varying unstable basic flow. It is seen from the GCM simulations, that the Himalayan orography plays a crucial role in anchoring the El Niño induced extratropical westerly troughs far to the west in the high latitude belt. In the absence of orography it is seen that the ENSO induced extra-tropical cyclonic anomalies tend to intrude southward into the monsoon region thereby destroying the regional scale circulations completely. Another effect due to the Himalayas is to generate lee waves on the eastern side of the topographic barrier which encircle the globe in the subtropics and midlatitudes.  相似文献   

Upper mantle structure of the Apulian plate from Rayleigh waves     

G. Calcagnile  G. F. Panza 《Pure and Applied Geophysics》1980,118(2):823-830

Phase velocities of Rayleigh waves for the Adriatic Sea area are obtained in the period range 25–190 sec along the path (l'Aquila-Trieste) AQU-TRI and 20–167 sec along the path (Trieste-Bari) TRI-BAI.The phase velocities are systematically higher than the known values for the surrounding regions. The data inversion indicates the presence of a lithosphere typical of stable continental areas with clear high-velocity lid (V _s 4.6 km/sec) overlying a well developed low velocity zone (V _s 4.2 km/sec).P. F. Geodinamica C.N.R., Roma Pubbl. N. 189.  相似文献   

A Lagrangian stochastic model for the trajectories of particle pairs and its application to the prediction of concentration variance within plant canopies     

A. M. Reynolds 《Boundary-Layer Meteorology》1998,88(3):467-478

A simple Lagrangian stochastic model for the trajectories of particle pairs in high Reynolds-number turbulent flows is presented. In this model, the velocities of particle pairs are initially correlated but subsequently each particle moves independently. The independent single-particle trajectories are simulated using Thomson's model [J. Fluid Mech. 180, 529–556, 1987]. This two-particle model exactly satisfies the well-mixed condition for Gaussian turbulence when length scales, characterizing the two-point Eulerian velocity correlation function, vanish. Temperature variances, due to heat released as a passive scalar from an elevated plane source, within a model plant canopy (Coppin et al. Boundary Layer Meteorol. 35, 167–191, 1986) are shown to be well predicted by the model. It is suggested that for strongly inhomogeneous flows, the two-point Eulerian velocity function is of secondary importance in determining the simulated trajectories of particle pairs compared to the importance of ensuring satisfaction of the two-to-one constraint (Borgas and Sawford. J. Fluid Mech. 279, 69–99, 1994); i.e ensuring that one-particle statistics obtained from the two-particle model are the same as those obtained from the corresponding one-particle model. Limitations of this modelling approach are discussed.  相似文献   

地震波在孔隙介质中低频衰减现象的粘弹性特征分析及近似（英文）     

凌云  韩立国  张益明 《应用地球物理》2014,11(4):355-363

由于介观尺度的孔隙流体流动,弹性波传播过孔隙岩层时在地震频段表现出较强的频散和衰减。Johnson理论给出了在任意孔隙形状的条件下,部分气水饱和孔隙介质的理论相速度和品质因子的解析解。本文在Johnson模型的基础上,通过对Q值曲线的低频和高频近似,推导了Q值曲线的近似公式,以及基于孔隙介质基本地球物理参数和孔隙斑块几何形态参数T和比表面积S/V的最大衰减Qmin近似公式。通过与理论值的对比,对Qmin近似公式存在的线性误差进行改正,进一步提高了精度。复杂的斑块形态对最大衰减Qmin和过渡频率ftr的都产生一定影响,且对ftr影响更大。因为数值模拟直接求解介观尺度的Biot孔隙介质方程需要极大的计算量,我们使用Zener模型建立了等效粘弹模型,有效地模拟了地震频带内的衰减和频散现象。  相似文献   

间断有限元方法的数值频散分析及其波场模拟          下载免费PDF全文

贺茜君  杨顶辉  吴昊 《地球物理学报》2014,57(3):906-917

数值求解波动方程是大尺度正演波场模拟、基于波动方程的地震偏移和反演成像的关键.本文针对求解二维声波方程的Runge-Kutta 间断有限元（RKDG）方法的数值频散问题,从理论推导和数值分析的角度进行了深入研究,并将其与近似解析离散化方法（Optimal Nearly Analytic Discrete Method,简称ONAD 方法）、Lax-Wendroff 修正方法、交错网格（Staggered-Grid,简称SG）方法的数值频散进行了比较研究.结果表明：RKDG方法以及近似解析离散化方法在压制数值频散方面要好于上述其他方法,特别是空间精度为3阶的RKDG方法,即使当空间步长取波长的一半,即一个波长内取2个网格点时,最大的频散误差也不超过1.67%.同时,我们也通过波场模拟对比研究了不同数值方法的数值频散问题,进一步直观地验证了数值频散的理论分析结果.  相似文献   

Reconstructing historic Glacial Lake Outburst Floods through numerical modelling and geomorphological assessment: Extreme events in the Himalaya          下载免费PDF全文

Matthew J. Westoby  Neil F. Glasser  Michael J. Hambrey  James Brasington  John M. Reynolds  Mohamed A. A. M. Hassan 《地球表面变化过程与地形》2014,39(12):1675-1692

Recession of high‐mountain glaciers in response to climatic change frequently results in the development of moraine‐dammed glacial lakes. Moraine dam failure is often accompanied by the release of large volumes of water and sediment, termed a Glacial Lake Outburst Flood (GLOF). Chukhung Glacier is a small (~3 km²) receding valley glacier in Mt. Everest (Sagarmatha) National Park, Nepal. Unlike many Himalayan glaciers, which possess a thick mantle of supraglacial debris, its surface is relatively clean. The glacier terminus has receded 1.3 km from its maximum Holocene position, and in doing so provided the space for an ice‐contact moraine‐dammed lake to develop. The lake had a maximum volume of 5.5 × 10⁵ m³ and drained as a result of breaching of the terminal moraine. An estimated 1.3 × 10⁵ m³ of material was removed from the terminal moraine during breach development. Numerical dam‐breach modelling, implemented within a Generalised Likelihood Uncertainty Estimation (GLUE) framework, was used to investigate a range of moraine‐dam failure scenarios. Reconstructed outflow peak discharges, including failure via overtopping and piping mechanisms, are in the range 146–2200 m³ s^‐1. Results from two‐dimensional hydrodynamic GLOF modelling indicate that maximum local flow depths may have exceeded 9 m, with maximum flow velocities exceeding 20 m s^‐1 within 700 m of the breach. The floodwaters mobilised a significant amount of material, sourced mostly from the expanding breach, forming a 300 m long and 100 m wide debris fan originating at the breach exit. moraine‐dam. These results also suggest that inundation of the entire floodplain may have been achieved within ten minutes of initial breach development, suggesting that debris fan development was rapid. We discuss the key glaciological and geomorphological factors that have determined the evolution of a hazardous moraine‐dammed lake complex and the subsequent generation of a GLOF and its geomorphological impact. © 2014 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.  相似文献